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Abstract 



Two-way relaying is a promising technique to improve network throughput. However, how 
to apply it to a wireless network remains an unresolved issue. Particularly, challenges lie 
in the joint design between the physical layer and the routing protocol. Applying an ex- 
isting routing protocol to a two-way relay network can easily compromise the advantages 
of two-way relaying. Considering routing path selection and two-way relaying together 
can be formulated as a network optimization problem, but it is usually NP-hard. In 
this paper, we take a different approach to study routing path selection for two-way 
relay networks. Instead of solving the joint optimization problem, we study the funda- 
mental characteristics of a routing path consisting of multihop two-way relaying nodes. 
Information theoretical analysis is carried out to derive bandwidth efficiency and energy 
efficiency of a routing path in a two-way relay network. Such analysis provides a frame- 
work of routing path selection by considering bandwidth efficiency, energy efficiency and 
latency subject to physical layer constraints such as the transmission rate, transmission 
power, path loss exponent, path length and the number of relays. This framework pro- 
vides insightful guidelines on routing protocol design of a two-way relay network. Our 
analytical framework and insights are illustrated by extensive numerical results. 



1 Introduction 



Two-way relay channel (TWRC) improve throughput by exploiting bi-directional inter- 
ference pp. A typical model for TWRC contains three nodes, as shown in Fig. HJ where 
A and B want to exchange data via relay Ri, assuming that all nodes operate in a half- 
duplex mode. Taking amplify-and-forward (AF) TWRC as an example, A and B transmit 
their packets simultaneously to relay in the first time slot. Then Rx amplifies and broad- 
casts the superimposed waveforms in the second time slot. After receiving it, A and B 
subtract its own signal to obtain their desired data. By utilizing interference instead of 
regarding it as noise, TWRC enables A and B to exchange data in two time slots, which 
is only half of the time needed in the conventional routing scheme. Therefore, two-way 
relaying has become a promising technique to improve network performance and has been 
extensively studied recently, e.g., pQ - [T2]. Although there are different two-way relay 
techniques pQ [5] [ID], this paper is focused on AF TWRC for the following two reasons. 
Firstly, AF is simple to implement and insensitive to the environment change, such as 
variations in the coding schemes or the modulation methods [5] [9]. Secondly, the perfor- 
mance of AF is comparable with other techniques. In [1], we see that at low transmission 
power, it has acceptable performance; at high transmission power, its performance is even 
better than most techniques. 

In spite of the advantages of TWRC, how to efficiently use this technique to acquire 
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Figure 1: Illustration of a TWRC. 

performance gains in a wireless network remains a challenging problem. For example, 
in Fig. [21 where A and B want to exchange data and there are more than one relays in 
between, if we want to transmit through AF TWRCs, one crucial question is that how to 
select an end-to-end routing path in the presence of TWRCs. To the best of our knowl- 
edge, no research work is reported on end-to-end routing in an AF two-way relay network. 
Thus, we are motivated to tackle this problem. In this paper, we develop a framework 
to analyze and compare the performance of different routes, from a perspective of infor- 
mation theory. Therefore, the three most fundamental parameters, i.e., energy efficiency 
(EE), bandwidth efficiency (BE), and latency, are considered as the metrics for evaluat- 
ing the performance of a routing path. Conventional routing metrics, such as hop count, 
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Figure 2: Illustration of a wireless network: A and B want to exchange data, but there 
are more than one relays in between. How to select an end-to-end routing path so that 
the route offers best performance in the presence of AF TWRCs? 
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ETX [13], ETT 1 14] . are not applicable to our framework. EE and BE are measures of 
how efficiently the network utilizes energy (including both transmission and processing 
energy) and bandwidth to transmit data. In this paper, we use a model for processing 
energy instead of assuming it as a constant as in [15J - [17J. Therefore, we are able to 
see with more accuracy how processing energy influences the network performance. Note 
that BE also characterizes the end-to-end rate in bit per channel use. Latency denotes 
the time elapse when data traverse the network. In other words, it measures the delay 
experienced by each bit before it reaches the destination. High EE, high BE, and low 
latency are the desired properties for a transmission scheme. However, normally they 
cannot be achieved simultaneously since tradeoff exists between them [18] [19]: we need 
to sacrifice one of them in order to gain the improvement in the other. 

The ultimate goal of this paper is to provide insights towards routing protocol design 
considering the fundamental requirements of BE, EE and latency. To build a framework 
for studying the performance of different routing paths from the perspective of BE and 
EE, we adopt a bottom-up method, beginning with analysis of a traditional three-node 
AF TWRC. Interestingly, given BE, we find a threshold for path loss exponent such that, 
when it is above the threshold, TWRC with relay located in the middle consumes the 
smallest energy. We also derive the power allocation associated with that condition of 
smallest energy consumption. After the simple case is studied, we extend the approach of 
analysis to an AF two-way relay network, where multiple pairs of nodes want to exchange 
data through more than one relays. With perfect scheduling, this two-way relay network is 
decomposed into a number of small networks each with only one source-destination (SD) 
pair. Thus we only need to consider one of those small networks, with an aim at selecting 
an optimal routing path for that particular SD pair. To better use the results gained in 
the previous analysis of three-node AF TWRC, we assume that each routing path has a 
small curvature and equi-spaced relays. Under this assumption, we formulate the three 
metrics, i.e., BE, EE, and latency, for routes with different number of relays. Numerical 
analysis is then carried out to analyze and compare the performance of different routes. 
Specifically, we find that the output parameters BE and EE of a route of a given SD pair 
are determined by the path length, number of relays and processing energy model. Since 
high EE, high BE, and low latency cannot be achieved simultaneously [18] [19], we define 
a general objective function integrating those three metrics. The optimal route can be 
found through this function. 

The following contributions are made in this paper: 

• We build an information theoretic framework to select an optimal routing path in 
an AF two-way relay network that provides the best tradeoff between EE, BE, and 
latency. 

• For AF TWRC, We find a threshold as a function of the transmission rate. When 
the path loss exponent surpasses the threshold, locating relay in the middle leads 
to the lowest energy consumption. 

• A power allocation scheme is developed to different routes. It allows each route to 
consume the smallest energy while still achieve the same transmission rate. 

This paper is organized as follows. Section II formulate the performance metrics 
for a traditional three-node AF TWRC, and analyzes how relay's location impacts the 
total energy consumption. Section III presents the network model and compares two 
schemes that enable multi-hop transmission through TWRCs. In Section IV we extend 
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the analysis frame to an AF two-way relay network, based on the Hop- by-Hop scheme. 
Section V studies the performance of different routes based on numerical results, and 
provides insightful observations towards routing protocol design. This paper concludes in 
Section VI. 

The discussions in this paper are based on the following notations. 

• Xf symbol transmitted by node i 

• yf symbol received by node i 

• Z{. additive white noise at node i, CAf(0,N ) 

• a: path loss exponent 

• /3: denote how many times relay amplifies the received signal 

• R: transmission rate at each link (bit/channel use) 

• Ri\ the ith relay 

• 7: bandwidth efficiency (bit/s/Hz) 

• R e : end-to-end rate (bit /channel use) 

• d route : length of a route (m) 

• Pf. transmission energy per channel use at node i (J/channel use) 

• Pproc,k'- average processing energy per channel use for a route with k relays (J/channel 
use) 

• £: energy efficiency (bit/ J) 

• k: total number of relays along a route, k = 0, 1, 6 

2 Three-Node Amplify-and-Forward TWRC 

In this section, we focus on a simple route with only one relay, i.e., a traditional three-node 
AF TWRC. We first introduce a model for computing processing energy, based on which 
the system's performance measures could be derived. We are then interested in finding 
when the route gives best performance, by varying relay's location and the amplification 
factor. Previous research has found that AF TWRC with relay in the middle gives the 
highest sum rate [3] and lowest outage probability [2D], with a given power allocation 
and fixed path loss exponent. But in this section we let both of the power allocation 
and path loss exponent be variables, and find that relay in the middle does not always 
provide the minimum energy consumption. One related work, i.e., [21], considers the 
power provision and relay placement problem for AF TWRC under Rayleigh fading and 
given outage probability, reaching a conclusion that relay is best positioned at the middle 
point to achieve lowest energy consumption. But our analysis is performed under different 
assumptions, i.e., large scale fading and capacity-achieving coding, which consequently 
leads to a different conclusion: whether relay in the middle is the most energy efficient 
depends on the relation between path loss exponent and transmission rate. 
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2.1 Model for Processing Energy 



When data is transmitted from source to destination, energy is dissipated at two places: 
in the radio hardware (i.e., processing energy, consumed at both source and destination) 
and in the EM waves (i.e., transmission energy, only consumed at the source node). For 
the processing energy, we assume a similar model as in [22], which is simpler than that 
in [22], yet still captures the main characteristics of the circuit power consumption. The 
total processing energy P proc for a transmitter-receiver pair is 

Transmitter: P T — ( l)Pt r + Pro, 

<7 (1) 

Receiver: P R = P R0 , 

where Pt and Pr denote processing energy dissipated at the transmitter and receiver; P tr 
is the transmission energy carried by the EM waves; rj is a constant representing the drain 
efficiency of the power amplifier (PA); P^ an d Pro are constants representing the energy 
consumption in radio electronics except PA, such as ADC, DAC, LNA, etc. 



2.2 Performance Measures 

For the derivation in this section and the following section, let Xi and i/i denote the 
transmitted and received symbols by node i, let be the noise at node i with a distribution 
of CAf(0, Nq), let Pi be the transmission energy per channel use at node i. Note that we 
do not distinguish between symbols transmitted or received in different time slots, for no 
ambiguity will occur during the derivation. Assuming that each node operates in half- 
duplex mode. As shown in Fig. [H A and B will transmit their packets simultaneously to 
relay in the first time slot,, i.e., 

y R = h A x A + h B x B + z R , (2) 

where h A and h B are the channel gains of the channels between A and relay, B and relay, 
respectively (assuming symmetric channels). After receiving the overlapped waveform, 
relay amplifies and broadcasts it in the next time slot, i.e., 

x R = (3(h A x A + h B x B + z R ), (3) 

where (3 denotes how many times relay amplifies its received signal. After receiving it, A 
and B will perform the so called "self-cancellation" to extract their desired data, i.e., 

Va ~ h A (3h A x A = h A (3h B x B + h A (3z R + z A , (4) 

y B - h B (3h B x B = h B f3h A x A + h B f3z R + z B . (5) 

Assuming that transmission is done by capacity-achieving codes and a common transmis- 
sion rate at each link. Then 

P = log 2 (l + 
= log 2 (l + 
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2.2.1 Latency 

A and B exchanges data every two time slots, so the latency experienced by each bit is 2 
time slots/bit. 

2.2.2 Bandwidth Efficiency (BE) 

Assuming that A and B exchange data with a common transmission rate R (bit / channel 
use) at each link, the end-to-end rate R e (bit/channel use) of the three- node system is 
R e = R. Given that each complex channel dimension occupies p seconds x hertz, where 
p is a constant depending on the techniques in the physical layer, then the BE (bit/s/Hz) 
of this system is 7 = R e / p [IB] . In this paper, we assume p — 1, which corresponds to the 
maximum BE 

7 = R e /p = R e = R. (8) 

2.2.3 Energy Efficiency (EE) 

Consider N total channel uses during the two time slots, then for AF TWRC, each time 
slot has N/2 channel uses. Let P proC) i be the average processing energy per channel use 
of the system, then its EE (bit/ J) is 

ReN 

^ ~ (P A + P B + P R )N/2 + P proC)l N- () 

According to ([1]), the total processing energy P proc ,iN can be divided into two parts: 
one is the processing energy in PA, which is linearly proportional to the transmission 
energy, the other is a constant value denoting energy dissipated in electronic circuits 
other than PA, i.e., 

PprociN = (I/77 - 1)(P A + P B + Pr)N/2 + P 0>1 N, (10) 

where (P A + Pb + Pr)N/2 is the transmission energy, and Po,i is the average circuit power 
consumption per channel use, except that consumed in PA. Substituting f[8l) ffT0]) into (J9]) 
gives 

_ 2R 
e ~ l(P A + P B + P R ) + 2P 0A 

2.3 Best Performance 

We are now interested in finding the highest EE of AF TWRC when BE is given. Seeing 
from (jSJ and (II ip . given BE, EE is maximized when the denominator in fill I) is minimized. 
Define / as the denominator of EE, i.e., 

f(R, h A , h B , (3) = l/r]{P A + P b + P R ) + 2P ,i. (12) 

/ equals the total energy consumption during two time slots, with each time slot occupying 
one channel use. Therefore, we call / the energy function. For simplicity, let hi = \h A \, 
h% = \hs\, and x = \(3\ 2 . Substituting (J6j) (JTJ) into f|T2|) gives: f(R,hi,h 2 ,x) = 

^,2(2H- 1 )(l^)( 1 + ^) + ^ +2poi (i3) 
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Since Po,i is a constant, and R is given by BE, there are three variables left, i.e., hi, h 2 , x. 
Thus, given BE, our goal is to find a vector (x*,h\,h 2 ) that offers the globally minimum 
value of /. 

Let d%, d 2 be the distance between A and relay, B and relay, respectively. For simplicity, 
we assume that the EM wave experiences large-scale fading only, i.e., h\ = h\ = d 2 a , 
where a is the path loss exponent. This is a common assumption when one wants to 
theoretically relate channel gain to distance (or location), e.g., [3] [15]. Assuming that 
the direct distance between A and B is fixed, denoted by d. If relay is put too far away, 
e.g., d\ > d, then A and B would prefer direct transmission without relaying. Therefore, 
the variables hi, hi and x are restricted by the following conditions 

x > 0, hi > 0, h 2 > 

hf /a < d, h 2 2/a < d. (14) 

h -y« + h f/« > d 

Taking partial derivative of / with respect to x, and let it be zero gives 



x 



2 < 2 ' i - 1 > (15) 



and df/dx < for x < xq, df/dx > for x > xq. Therefore, the global minimum point 
(x*, hi, h 2 ) of /, if exists, must satisfy (TT5]) . for if not, we can always find a Xq calculated 
from (fT5l) so that f{xQ,h\,h^) < f(x* , h\, h^), which contradicts the assumption that / 
achieves global minimum at (x*, h\, h 2 ). Substituting ffl5l) into f|T3l) gives 



ahih 2 + b(h\ + h\ 

J y-n, < b ii "-2, xq) = 

where 



f(R, h u h 2 , x Q ) = , 2 V , 2 1 lJ - + 2P ,i, (16) 

r]n 1 n 2 



a = 2N ^(2 R + 1 - l)(2 fl + 1 - 2) (17) 

and 

b = N (2 R+1 - 2). (18) 

Therefore, finding the minimum energy consumption is equal to finding the global 
minimum value of (TTB|) . Since df /dhi < and df/dh 2 < 0, given hi (h 2 ), f decreases 
when h 2 (hi) is increased. Then the minimization will be achieved when hi, h 2 (h 1 2 ^ a , 
h 2 2 ^ a ) are as large (small) as possible, so from (114)) . we have that / achieves its minimum 
value when hi 2 ^ a + h 2 2 ^ a = d. This means that TWRC consumes less energy when relay 
is located along the straight line connecting A and B, which agrees with our intuition. 
Let h l 2,/a = dcos 2 8, h 2 2 ^ a = dsin 2 8, where 9 e (0,7r/2), then ([TBI) becomes 

ad a cos a 9 sin a 9 + bd a (sin 2a 9 + cos 2a 6) 
J (9) — h 2P ,i- (19) 

Its first derivative is zero when 9 = 7r/4, which corresponds to the situation when 
relay is in the middle between A and B. To see whether / achieves local minimum at this 
critical value, we check its second derivative and get 



{local min a > 1 + „ „ , , 
local max a < 1 + 
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Transmission Rate at A or B (bit/channel use) 

Figure 3: The (R, a) region where energy consumption is minimized when relay is in the 
middle. 



From f )20|) . we know that whether relay in the middle achieves local minimum energy 
consumption is determined by whether the path loss exponent is above the threshold 
formed by the transmission rate. Fig. [3] shows in what (R, a) region relay in the middle 
consumes locally minimum energy, from which we see that this proposition holds true 
in most cases. In 802.11a standards [23], for example, we have 54 Mbps and 20 MHz, 
so R = 2.7 bit/channel use, then as long as a > 2.04, transmission through TWRC 
with relay in the middle saves the most energy. Fig. |4] plots the energy function divided 
by iVn for both AF TWRC and direct transmission between A and B, assuming R = 1 
bit/symbol, a = 2 - 2.4, d = 20 m, and P pr0 c,i — 0. The x-axis denotes a relay's location 
along the line connecting A and B. Note that for different d and P proc ^, the absolute value 
of each curve will be different, but here we are interested in their shape as well as the 
comparison results of their magnitudes, which will not change with d and P pr0 c,i- In ( 120]) . 
let R — 1 bit/symbol, we get that when a > 2.22, / is locally minimized when relay 
is in the middle, which can be confirmed by Fig. HI As shown in this figure, / achieves 
not only local minimum but also global minimum under the condition a > 2.22. When 
a < 2.22, / is locally maximized when relay is in the middle, although not always globally 
maximized, as shown in Fig. HI As a decreases to 2, the energy consumption for TWRC 
tends to be higher than direct transmission, which means that direct transmission works 
better in terms of less energy dissipation. Together with ( 12"U|) . we conclude that as long 
as a > 1 + v/ fSrEf ) resources are better utilized when transmitting through TWRCs, 
rather than direct transmission; and relay is best located at the middle point. Discussions 
in the following paper will be focused on this (R,a) region. 



3 AF TWRC in A Wireless Network 

Previously, we studied a simple route with only one relay, and derived its optimal EE 
given BE. In this section, we consider routes with more than one relays. We first propose 
the network model. To analyze different routes in terms of EE, BE and latency, we need a 
scheduling scheme that enables transmission through AF TWRCs along a multihop route, 
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Figure 4: Total energy consumption (normalized by Nq) for various locations of the relay, 
when R = 1 bit/symbol, a varies from 2 to 2.4, d = 20 m, r\ = 1, Pq.i = 0. Energy for 
direct transmission is plotted in straight lines. 

so a simple scheme is presented. Based on this simple scheduling scheme, we then derive 
the performance measures of routes with different number of relays. 

3.1 Network Model 











..... 





Figure 5: A and B want to exchange data through AF TWRCs, which route performs 
best? 

A single-frequency AF TWRC network contains multiple pairs of nodes that want to 
exchange data through multiple relays. Perfect scheduling is assumed in the network so 
that each source-destination (SD) pair transmit in its own time slot with no concurrent 
transmission inside its interference range. Therefore, the whole network can be decom- 
posed into multiple small networks each with only one SD pair. Moreover, considering the 
whole network is equivalent to considering each of the small networks separately. Thus, 
it's enough to analyze the routing path selection problem for one small network, i.e., one 
SD pair. Since our analysis is based on the assumption of perfect scheduling, the results 
gained will provide upper bound of the network performance. 

We now focus on one SD pair. Suppose A and B want to exchange data through 
TWRCs, as shown in Fig. |5j Among the several routes, each with different length and 
number of relays, we want to know which route gives the best performance. To better 
use the results in Section II, we make the following two assumptions: 
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• the relays along each route are equidistant; 

• each route has a small curvature so that it can be treated as an elongated straight 
line. 

The two assumptions are reasonable. Practically, the relays will not be allocated too 
near or too far from each other. Besides, based on the analysis from previous section, in 
most cases, relay located in the middle consumes the smallest energy, and hence equally 
spaced relays along a route tends to achieve low energy consumption. Route with a large 
curvature means that most energy is spent for data in moving around the source node 
rather than forwarding to its destination. Besides, nodes with several hops in between 
may be close in distance so that they cannot transmit simultaneously, which decreases the 
concurrent number of transmission nodes. Therefore, route with a large curvature needs 
to be avoided. 

In addition, let k be the number of relays along a route, this paper considers the case 
when < k < 6, for TWRC is usually used in a small-scaled wireless access network, 
which is sure to have an upper limit on the number of hops. Note that our framework can 
be extended to the case when k > 7, with a possibly increasing complexity as k becomes 
larger. For simplicity, we also assume that each node works in half-duplex mode and that 
each hop transmits at a constant rate. This constant rate control mechanism ensures 
that each node forwards a packet to the next hop per time slot, which help maintain the 
stability of the system. 

3.2 Two Multi-hop Transmission Schemes 

To enable multi-hop transmission through TWRCs, two schemes can be used: Hop-by- 
Hop scheme and End-to-End scheme. We assume that the two end nodes always have 
packets to exchange. Taking k = 5 as an example, Fig. [6] illustrates how Hop-by-Hop 
scheme works. At the beginning, no packet is traversing along the route. Then the end 
nodes start to put in the packets that they want to exchange with each other. After two 
time slots, the system enters a stable state in which packets are transmitted in a recursive 
pattern: during every 4 time slots, the end nodes will insert one new pair of packets into 
the system (i.e., packets that will be delivered to the other side), and receive one new 
pair of packets from the other side (i.e., packets that were inserted by the other side at 
an earlier time slot); all the relays in between will help forward data through TWRCs. In 
that way, the total time slots needed to exchange n pairs of packets is 2 + 4n, when k = 5. 
However, the recursive pattern is different when k is different. As is illustrated in the 
next section, when k is odd, then all the nodes are involved in TWRCs; when k is even, 
then one of the end node will perform unicast transmission. Compared with Hop-by-Hop 
scheme, which is a loose (4 time slots / exchange) pattern formed by concatenated three- 
node subsystems and unicast channels, End-to-End scheme is a compact (2 time slots 
/ exchange) pattern formed by overlapped three-node subsystems, where each node will 
interchange between transmitting and receiving every time slot. 

It seems that End-to-End scheme is better since it is simpler and has higher through- 
put, however, this scheme experiences severe noise accumulation between adjacent three- 
node subsystems. As mentioned previously, End-to-End scheme has a compact transmis- 
sion pattern, where only the two end nodes decode the received packets. Each relay in 
between performs subtraction, amplification, and broadcast. Correspondingly, the noisy 
packets will be subtracted from or superimposed with other noisy packets, resulting in 
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Figure 6: Illustration of the Hop-by-Hop scheme for k = 5. 



a rapidly growing noise as it traverses the network. Fig. [7] shows how noise accumulates 
at the end nodes as time elapses, assuming that the initial noise variance is 1, relay does 
not do amplification, and noise accumulation due to substraction is ignored. The expo- 
nentially increased noise will terribly distort signals received at the end nodes. As to the 
Hop-by-Hop scheme, noise will not be accumulated between adjacent subsystems, since 
the two end nodes in each subsystem will decode the packet each time the relay broad- 
casts the overlapped waveforms. Therefore, although Hop-by-Hop scheme spends twice 
the time for the end nodes to exchange one packet, its performance is still better than 
End-to-End scheme. 




Figure 7: In End-to-End scheme, noise received at the end nodes increases exponentially 
with the number of received packets. 



4 Performance of Hop-by-Hop Scheme 

In this section, we investigate the performance, i.e., BE, EE and latency for routes with 
different number of relays, based on the Hop-by-Hop scheme. Since how power is provi- 
sioned at each node affects the system's EE and BE, we will determine the optimal power 
allocation scheme under which a route achieves the highest EE with a given BE. The 
difficulty during derivation mainly comes from the interference influence. Since we are 
utilizing Shannon's capacity formula to derive the fundamental relation, any interference, 
no matter how small it is, should be included in the SINR part inside the formula, which 
increases complexity during the derivation process. Numerical analysis in the Section 
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Figure 8: Illustration of the Hop-by-Hop scheme: k = 0. 



IV illustrates that ignoring this interference will cause big error percentage at high BE. 
Therefore, we need to carefully tackle with this interference influence. 

Recall that relays are equally spaced, and that each node is transmitting at rate R 
(bit/channel use) with capacity-achieving codes. Besides, each node operates in half- 
duplex mode. Denote by h and d the channel gain and the distance between two consec- 
utive nodes (for simplicity, we use h as both the channel gain and its Euclidean norm). 
Here, we only consider large-scale fading, i.e., h 2 = d~ a . Let Pi denote the transmis- 
sion power per channel use at node i. We do not distinguish between h and Pi for 
different k, since there is no ambiguity in the following derivation. Let P pr0 c,k be the 
average processing energy per channel use for route with k relays, then Section II gives 
P P roc,k = (1/v — ^)Ptr,k + Po,k- Assuming that each channel use occupies 1 second x hertz. 
Besides, let transmission range be 1 hop, interference range be 2 hops, i.e., nodes within 
two hops cannot transmit simultaneously, except in TWRCs. 



4.1 Latency 

In Hop-by-Hop scheme, packet is forwarded to the next hop per time slot, so the latency 
of a route with k relays is k + 1 time slot/bit. 



4.2 k=0 

Fig. [8] illustrates k = 0, i.e., direct transmission. 

• BE and EE 

2R 

7 = jR ' e = u(R,h) + 2P procfi (21) 

u(R, h) = Pa + Pb is the transmission energy during the two time slots, with each time 
slot occupying one channel use. Besides, 2P proCt0 = (1/rj — l)u(R, h) + 2P ,o- 

• Optimal Power Allocation 

Given BE, u(R, h) achieves its minimal value when using capacity-achieving codes, 
i.e., ^ 

i2 = log 2 (l + ^*), z = A or B. (22) 

iVO 

Thus, the optimal power is 

Pa = Pb = N (2 r - l)h- 2 . (23) 
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4.3 k=l 



Fig. [T] depicts the case when k = 1, i.e., three-node TWRC. As analyzed in Section II, A 
and B exchange one packet every two time slots. 

• BE and EE 

2R 

^ ' ^ h) + ZPproc,! ^ 

g(R, h) = Pa + Pr 1 + Pb is the transmission energy during the two time slots, with each 
time slot occupying one channel use. 1P vroC) \ = (^/v ~ 1)#C^j h) + 2Po,i- 

• Optimal Power Allocation 

Given BE, the minimal value of g(R,h) is g(R,h) m i n = (a + 2b)h~ 2 , which can be 
derived from (FTO]) by setting hi = h 2 = h, 77 = 1, and Po,i — 0, with a and b given in 
(TT7I) (TT5|) . The corresponding optimal power allocation can be found by setting tiA = h B = 
h in ©©, i.e., 

Pa = Pb = (2 r - l)(h 2 p 2 + l)N /(h A p 2 ), 

P R = (2h 2 P A + N )(3 2 . [Zb) 

The optimal f3 2 is given in (TT5]) . 



4.4 k=2 

Fig. M illustrates how Hop-by- Hop scheme works when k = 2. During the first and fourth 
time slots, A sends and receives one packet from R4 through unicast transmission. During 
the second and third time slots, B sends and receives one packet from R 2 through TWRC. 
In this way, they form a recursive pattern of 4 time slot /exchange. 

• BE and EE 

2R 

7 = jRA ^ = u(R,h)+g(R,h) + 4P proc , 2 (26) 
u(R, h) = Pa + Pr x and g(R, h) = Pr 1 + Pr 2 + Pb represent transmission energy through 
unicast channel and TWRC, respectively. 4P pr . OCj2 = (1/rj — l)(u(R, h) + g(R, h)) +4P 0j2 . 

• Optimal Power Allocation 

Given EE, the minimum of u(R, h) and g(R, h) have been found in previous cases of 
k = and k — 1, with the optimal power provision given by fl23l) fl25l) . 

4.5 k=3 

Fig. QUI presents the Hop-by- Hop scheme for k = 3, where the left half nodes and the right 
half nodes form two TWRCs, respectively. 

• BE and EE 

2R 

7 = jR/2 ' e= gi(R,h)+g 2 (R,h)+4P proc , 3 (2?) 
gi(R,h) = Pa + Pr 1 + Pr 2 and g 2 (R,h) = Pr 2 + Pr 3 + Pb represent the transmission 
energy for the two TWRCs. 4P proCi3 = (I/77 - l)(#i(P, h) + g 2 (R, h) + 4P , 3 . 

• Optimal Power Allocation 

Given EE, the minimal values of gi(R,h) and g 2 (R,h) are the same, as given in 
previous case of k — 1, with the optimal power allocation in <^^. 
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Figure 10: Illustration of the Hop-by-Hop scheme: k 



3. 



4.6 k=4 

Fig. HH depicts the case for k — 4. During the first two time slots, data is exchanged 
between R2 and R4. In the third time slot, A inserts a new packet into the system while 
B receives a packet. After that, the role of A and B exchanges in the fourth time slot. 
Thus, a recursive pattern is built with 4 time slot /exchange. 



BE and EE 



2R 



(2* 



g(R, h)+t(R, h)+AP procA 

g(R, h) = Pr 2 + Pr 3 + Pr 4 , t(R, h) = Pa + Pr 1 + Pr 2 + Pr 4 + Pb are the transmission 
energy consumed during the first two and last two time slots, respectively. 4P procA = 
(l/ V -l)(g(R,h) + t(R,h))+4P 0A . 

• Optimal Power Allocation 

Given EE, the minimal value of g(R, h) have been derived in the case of k — 1, where 
Pr 2 i Prz an d Pr 4 are given in f[25l) . Now we will find the minimal value of t(R, h). A, Ri, 
R 2 form a TWRC, i.e., 



VRi 



hx A + hx R2 + z Rl + V3~*hx Ri . 



hxR 1 + za + Vo a hxi 



VA = 

Vr 2 = hx Rl + zr 2 + V3~*hx B . 
B, R 4 perform unicast transmission, i.e., 



Vb 



hx R , + z B + V3 a hxR 2 + V5 



y hx A , 



(29) 

(30) 
(31) 
(32) 

(33) 



Figure 11: Illustration of the Hop-by-Hop scheme: k = 4. 
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y Ri = hx B + z Ri + V3^hx Rl . (34) 
If capacity-achieving codes are used, then 



3- a h 4 fi 2 P R4 + {h 2 (5 2 + 1)N + 5~ a h 2 P B 



(35) 



-2 fl -l, (36) 



3- a h 4 /3 2 P R4 + (h 2 /3 2 + 1)N + 3- a h 2 P B 



1,2 p 

11 ^ R * = o R - 1 (371 

3- a h 2 P R2 +5- a h 2 P A + N ' 1 ' 

h2PB -2 R -1. (38) 



From (1301) we also have 

P Rl = (3 2 (h 2 P A + h 2 P R2 + 3- a h 2 P Ri + N ). (39) 

There are 5 equations, i.e., (|35l) - (l39j) . and 5 variables, i.e., Pa, P Rl , Pr 2 , P Ri , P B - Thus, 
we are able to solve the equations and get the 5 variables. 
(EHDdSnD gives 

P Rl = - 3^ 2 iV . (40) 

Substituting (JIT?]) and (13T1) into the formula of t(R, h) gives 

t = (l + (2 R -l)5- Q )P A +(l+(2 R -l)3- a )P i?2 + (^ T + l)P B +((2 i? -l)-3 Q )/i- 2 iV . (41) 

Eliminating P Ri from fl35|) and f )38|) gives 

2 R h 4 p 2 P R2 + h i f3 2 {2 R - l)P A = (3 a h 2 + 5- a h 2 (2 R - l))P B + (2 R - 1)(1 - 3")AT . (42) 
Eliminating Pr 4 from (136]) and (1381) gives 

2 R h i f3 2 P A + h i (3 2 (2 R - 1)P* 2 = (3 a /i 2 + 3- a h 2 (2 R - 1))P B + (2 R - 1)(1 - 3 a )N . (43) 

From (USD and (g3]), we get 

(2 R+1 - l)h 4 (3 2 P A = Bh 2 P B + (2* - 1)(1 - 3 a )AT , (44) 

(2 fi+1 - l)/i 4 /3 2 P R2 = a 2 h 2 P B + (2 R - 1)(1 - 3 Q )iV . (45) 
Substituting (J37J) into ([38]), together with (jUJ) and (T45]), we get 

(1 - a 5 )h 2 P B = (o 3 + /* 2 /3 2 a 4 )iV . (46) 
Substituting (pH]) (135D (|35j l into the formula of (JH]) gives 

t(P, /i) = (a 8 /3 2 + ^ + a e h- 2 )N > (2^ES + ae)^^" 2 . (47) 

So the minimum of t is t(R, h) min = (2^ya 7 a s + a 6 )A^ /i~ 2 , achieved when (3 A = a 7 h~ 4 /a 8 . 
The corresponding optimal power allocation is 

(q 3 + a 4 fr 2 /3 2 )iVo 
S = (l-a 5 )tf ' (48) 
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Figure 12: Illustration of the Hop-by-Hop scheme: k = 5. 



ai h 2 P B + (2 R -l)(l-?> a )Nv 



R2 



(2 R + X - l)/i 4 /3 2 
a 2 h 2 P B + (2 R -l)(l-3 a )iV 



(49) 



(50) 



(2*+* - l)/i 4 /3 2 

The optimal power allocated at Ri and R4 can then be easily computed from f )37|) fl38 
{aj}f =1 are functions of R and a: 



a x = -{2 H - 1) 2 5~ Q + 2 K (2 H - 1)3~ Q + 3 Q , 
a 2 = 2 R (2 R - l)5- a - (2 R - l) 2 3- a + 3 a , 
a* = ^ (1 - 3 Q )(3-°2 + (45~ a + 27~ Q )(2 fi - 1)) + (2 R - 1) 



a 5 



2 R+i - 1 



2 R+i - 1 



a 4 = 9~ a {2 R - 1) 2 + 3~ Q (2 K - 1) 



a/ n R 



a 2 (3-° + 45- Q (2 fi - 1)) + ai(3~ Q + 27~ a (2" - 1))] 



a/ n R 



a 6 



a 7 



a 2 a 4 (l + (2 R - 1)5-) + o ia4 (l + (2 fl - 1)3'°) | fl ga | a 3 



3° 



(2«+i_i)(i_ a5 ) 



l-a 5 K 2 R -l 



(51) 
(52) 
(53) 
(54) 
(55) 

+ 1), 
(56) 



(2 R - 1)(1 - 3°) 
2«+! - 1 



(2+(2 R -l)(3- Q +5- a )H 



a 4 , 3 a 



a 2 a 3 (l + (2* - 1)5-") a^jl + (2 K - l)3~ a ) 

(l-a 5 )(2*+i-l) + (l-a 5 )(2«+i-l) 

(57) 



1 - a 5 K 2 R - 1 



!)■ 



(5? 



4.7 k=5 



Fig. [12] illustrates how Hop-by-Hop scheme works when k = 5. This recursive pattern has 
been analyzed in previous subsection. 



BE and EE 



2R 



(59) 



g(R, h) + s(R, h) + 4P proc , 5 

g(R, h) = Pr 2 +P R3 + Pr 4 , s(R, h) = P A + P Rl +Pr 2 + Pr 4 +Pr b +Pb are the transmission 
energy consumed during the first two and last two time slots, respectively. 4P proCi5 = 
(l/rj-l)(g(R,h) + s(R,h))+AP 0j5 . 



Optimal Power Allocation 
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Similar to the previous case, given EE, the minimal value of g(R, h) has been derived 
when k = 1. Next we find the minimal value of s(R,h) and the corresponding power 
allocation. A, R 1; R 2 , and B, R 4 , R 5 form two TWRCs, i.e., 

y Rl = hx A + hx R2 + V3~*hx R4 + \f^r° t hx B + z Rl , (60) 

x Rl = fiy Rl , (61) 

x A = hx Rl + V5-^hx Rs + z A , (62) 

xr 2 = hx Rl + V¥^hx R5 + z R2 . (63) 
Assuming capacity-achieving codes, then 



R2 _ nR 



h 2 (P Rl - Ph\P A + P R2 )) + h-«WP Rb + N 



h\P Rl - f3 2 h 2 (P A + P R . 2 )) + 3- a h*P Rs + N 
From flBTl). we also have 



1, (64) 
1. (65) 



P Rl = P 2 (h 2 P A + h 2 P R2 + 2r a h i P IU + 5~ a h 2 P B + N ). (66) 

Due to symmetry, P A = P B , Pr 2 = P Ri , P Rl = P Rs , so we only need to consider the one 
of the two TWRCs. Accordingly, s(R,h) = 2{P A + P Rl + Pr 2 ). Similar to the case of 
k = 4, we can solve the 3 equations, i.e., (164]) - (16"6"|) . 
Substituting (I6"6]) into the formula of s(R, h) gives 

s = 2(1 + (1 + 5~ a )h 2 (3 2 )P A + 2(1 + (1 + 3- a )h 2 f3 2 )P R , 2 + 2(3 2 N . (67) 

(16151) minus ([54]) gives 

h A p 2 (P A - P R2 ) = (2 R - 1)(3~ Q - 5- a )h 2 P Rl . (68) 

Substituting (I6"B"|) into (I6"%j) gives 

h 4 f3 2 P A = h^^PR, + b 2 h 2 (5 2 N . (69) 

Substituting (169]) into (164]) and ([6"o]) gives 

h A (3 2 P R2 = (b 3 h 2 (3 2 + bt)N . (70) 

h^ 2 P A = ({hb 3 + b 2 )h 2 (5 2 + 6!& 4 )iVo. (71) 
Substituting ([TO]) and (T7T1) into (1671) gives 

/i) = (6 5 /3 2 + ^ + b 7 h- 2 )No > (2^b 6 + & 7 )Ar ^ 2 . (72) 

So the minimum of s is s(R, h) m m = {2^/b^ + b 7 )N h~ 2 , achieved when /3 4 = b%h~ A jb^. 
The corresponding optimal power allocation is 

p ((M 3 + fr 2 )fr 2 /3 2 + WiVo 

Pa = ^ ' ( 73 ) 
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p* = 13 uir ° (74) 



(hh 2 (3 2 + 6 4 )iV 

The optimal power allocated at Ri can be easily computed from (166)) . {bi}J =1 are compli 
cated functions of R and a, and their full expressions are 



l + (2 f? -l)(3- Q -5-")(l + 3-") 
1 l-(2«-l)(3-°-5- a )(l + 5-«)' 1 J 



(2 fi -l)(3- a -5- Q ) 



1 - (2 R - 1)(3- Q - 5-«)(l + 5- Q ) 



(76) 



, [(5-°2 + 25- Q )& 2 + 5~ Q + 1](2* - 1) 

3 l-(2 /? -l)[6 1 5- Q (2 + 5- Q ) + 5- a + 3- a + 15- a ]' 1 j 

2* - 1 

64 ~ 1-(2 R - I)[6i5- Q (2 + 5~ Q ) + 5~ a + 3~ a + 15~ a ] ' 

65 = 2(6x63 + 6 2 )(1 + 5- a ) + 26 3 (1 + 3- a ) + 2 (79) 
6 6 = 2(6 X + 1)64 (80) 
6 7 = 2(6x6 3 + 6a) + 26164(1 + 5~ Q ) + 26 3 + 26 4 (1 + 3- Q ). (81) 



4.8 k=6 



The recursive pattern of k = 6 resembles that of k = 4 and k = 5. During the first two 
time slots, the system works like in the third and fourth time slots when k — 4. During 
the last two time slots, the recursive pattern is the same as that when k = 5. 

• BE and EE 

9 R 

t(R, h) + s(R, h) + 4P pr . OCj6 

t(R, h) = P R2 + P R3 + P Ra + P R6 + P B and s(R, h) = P A + P Rl + P R , 2 + P Ri + P Rb + P Ra 
give the transmission energy during the first and the last two time slots, respectively. 
4P proc , 6 = (l/v - l)(t(R, h) + s(R, h)) + 4P ,6- 

• Optimal Power Allocation 

Due to the same recursive pattern, the minimal values of t(R, h) and s(R, h) as well 
as the optimal power allocation have been derived in previous cases of k — 4 and k = 5. 



5 Numerical Results 

In Section VI, we have derived the performance measures of routes with different number 
of relays. We also found the optimal power allocation scheme, under which the route 
consumes the smallest energy while still transmits at the same end-to-end rate. In this 
section, we numerically study the performance of different routes with the application of 
the optimal power allocation scheme. 

Let path loss exponent a = 4, then from (l20i) we know that only when the transmis- 
sion rate R > 0.087 bits/channel use, does the relay in the middle minimize the whole 
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Figure 13: Comparison of the EE-BE relation for direct transmission (k = 0) and TWRC 
(k = 1) when end-to-end distance d = 1000 m, a = 4, and processing energy is ignored. 



energy dissipation. Besides, set noise variance as N = —174 dBm/Hz, and drain effi- 
ciency f] = 0.75, which is an achievable value for high-class power amplifiers |22j. For the 
processing energy, we assume that Po,fc is proportional to the average number of senders 
and receivers per channel use, i.e., Po,k/Po,o — rn/2, where m is the average number of 
active nodes per channel use for route with k relays. Let Po,o = 5 x 10~ 7 mJ/channel 
use, then P ,i = 3/2P , , P),2 = 5/4P , , P , 3 = 7/4P , , P , 4 = 2P , , P , 5 = 9/4P , , 
P 0i 6 = 11/4P 0)0 . d route is the total length of the route. 

Fig. [13] compares the performance for direct transmission and TWRC, which corre- 
sponds to k = and k = 1. Since both of them spend 2 time slots in exchanging one pair 
of packets, they have the same latency, so only EE and BE are compared in the figure. 
Here we ignore the processing energy, since nonzero processing energy will not influence 
the results. For direct transmission, EE and BE are always negatively related, which is a 
fundamental feature inherited from the Shannons capacity formula [T3]. However, TWRC 
changes the way that EE-BE interacts: EE and BE are positively related at low rate; and 
then involved in a tradeoff relation after EE reaches its maximum value. In other words, 
at low rate, we can decrease energy consumption and increase transmission rate at the 
same time. This feature can also be seen from the equation (|24p . If processing energy is 
ignored, then EE for a three-node TWRC is 

£ = / = • (83) 

N {J{2 R + 1 - l)(2 R+1 -2) + 2 R + 1 - 2)h- 2 



It can be verified that lim^o £ — and lim^oo £ = 0, there must be a transmission 
rate that achieves maximum EE. From Fig. [131 we can estimate that this rate is R ~ 0.6 
bit/channel use. Besides, from the figure, we see that when rate approaches zero, direct 
transmission tends to have higher EE than TWRC, which is consistent with analysis in 
Section II. Therefore, in the following analysis, we focus on R > 0.087 bit/channel use, 
where TWRC will have better performance than direct transmission, and energy con- 
sumption is minimized when relay is in the middle. 



In Fig. 14(a) , we plot EE and BE along a route with length 1000 m and k varies 



from to 6. Again, we ignore the processing energy here, so only transmission energy is 
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Figure 14: EE-BE relation for different number of relays when a = 4, processing energy 
is ignored, (a) d route = 1000 m, k varies from to 6 (b) d route = 1000 m or 1200 m, k 
varies from 3 to 6. 



considered in this figure. Given the length of a transmission path and a certain BE, EE 
will increase with increased number of equidistant relays. In other words, if two routes 
connecting A and B have the same length, then the route with more relays consumes less 
transmission energy while still achieves the desired end-to-end rate (i.e., BE). This can 
be explained as follows. Given the length of a path, more relays along that path means 
smaller distance in each hop, and hence smaller transmission energy consumption in each 
hop. This effect dominates the increase in energy consumption resulting from more relays 
are consuming energy. 

Fig. 14(b) illustrates how length of the route affects its performance. Given BE and 
the total number of relays, EE decreases as the length of the route increases, since longer 
distance needs larger transmission power to maintain the same transmission rate. Combin- 
ing Fig. 14(b) and Fig. 14(a) , we conclude a general principle when comparing different 



routes in terms of EE and BE: given end-to-end transmission rate, route with shorter 
length and more relays tends to consume less energy. This principle has a key assumption 
that processing energy is small compared with transmission energy. But this is not always 
true: processing energy may dominate the total energy consumption at low transmission 
rate. 

Fig. 15(a) and Fig. 15(b) depict EE and BE relation for P ,o = 5 x 10~ 6 and Po,o = 
5 x 10 -7 mJ/channel use, respectively. With a high processing energy, e.g., in Fig. 15(a 



given the length of the route, a route with two relays in between has the best EE for most 
transmission rates. The main reason is that P pr0 c,2 is smallest among different P pr0 c,kS, 
k = 1,...,6. If processing energy is small, e.g., in Fig. 15(b), then the EE-BE relation 



resembles that in Fig. 14(a) , in the sense that a route with more relays tend to has larger 
EE as BE gets higher. Besides, if the number of the relays are given, the transmission 
rate which gives the maximal EE is different for different k. When k = 2, the system can 
achieve best EE by transmitting at an end-to-end rate around 1.1 bit/channel use. 



Fig. dH explains the reason why we cannot ignore the interference influence when using 
capacity formula to derive EE and BE relation. The error percentage when using SNR 
instead of SINR to compute capacity increases from 4.1% at BE= 1.2 bit/s/Hz to 325% 
at BE= 2.75 bit/s/Hz. Except the large difference in the numerical results of EE, there 
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Figure 16: Comparison of EE-BE relation with (SINR) and without (SNR) considering 
interference: d route = 1000 m, a — 4, k — 4, processing energy is ignored . 



20 




2 3 4 

Total Number of Relays 



Figure 17: Upper limit of the transmission rate and latency for various number of relays 
where d route = 1000 m, a = 4. 



is another drawback when we only consider SNR, that is, it cannot present the upper 
limit of the achievable transmission rate. If we use SNR in the capacity formula, we can 
achieve any rate as long as the transmission power is large enough, however, in reality it 
is not true. Taking k = 4 as an example, in (|35|) . Pr 2 , Pr 4 and Pb are of the same order, 
so there is an upper limit of R on the left-hand side. 

In Fig. [T7] we plot the upper bound of the transmission rate as well as latency for 
different ks. The reason for bounded transmission rate is that the interference from other 
nodes is of the same order as the transmission power, resulting from the common rate 
transmission scheme. Since there is no interference in the transmission pattern when 
k < 3, no upper limit of rate exists. Previous analysis, i.e., Fig. 14(a), shows that 



more relays along a path can decrease the total transmission energy, but there is one 
big disadvantage associated with more relays, that is, latency. For a route with k relays, 
latency is k + 1 time slot/bit, since each bit will be forwarded by k + 1 hops before 
it reaches destination. It is different from the time which the two end nodes spend in 
receiving consecutive packets from each other, which is 2 time slots for k < 1 and 4 time 
slots for k > 2. 



Previous analysis, i.e., Fig. 14(a) , tells us that more relays of a path tends to decrease 
the total transmission energy. However, more relays incurs larger latency. Thus, we need 
to jointly consider the three performance metrics, in order to decide which one of the 
three routes in Fig. [5] performs best. A route with high EE, high BE and low latency is 
preferred, so a general performance metric can be defined as 



EE/EE max x BE/BE max = EE x BE 
latency/latency max latency 



x K, 



where EE max , BE max , and latency max are the maximum achievable EE, BE, and latency of 
a system, e.g., Fig. |5] and K is a constant formed by the three maximum values. F is an 
objective function that we want it to be as large as possible. Note that the definition of 
F is not unique, e.g., if one cares more about EE, a different F function can be designed 
so that route with a large EE tends to have a high F value. 



Fig. 18(a) and Fig. 18(b) depict F/K for the three routes in Fig. |5j assuming their 
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Figure 18: Comparison of the three routes in Fig. [5j a 
(b) P ,o = 5 x 10- 7 . 



(b) 



0.75, 



Pi 
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configurations (length and number of relays) are 1200 m with k — 1, 1000 m with k = 
3, and 1600 m with k = 2. When processing energy plays a dominant role, e.g., in 
Fig. 18(a), the first route has the best performance, since it has the lowest latency and 



small processing energy consumption. When processing energy is smaller compared with 



transmission energy, e.g., in Fig. 18(b), the advantage of the second route becomes more 
obvious as the transmission rate increases. If processing energy can be ignored, according 
to the general principle weve analyzed previously, the second route has the best BE-EE 
relation since it has both the shortest length and the largest number of relays. 



6 Insights On the Routing Path Selection 

6.1 Summary of the Previous Results 

• Given BE, whether relay in the middle point achieves the highest EE is decided by 
(|2"01 . When the path loss exponent is larger than the threshold, relay is best located 
in the middle, besides, TWRC is more energy efficient than direct transmission. 

• TWRC, unlike normal channels, displays a different EE and BE relation: they are 
positive related at low BE and then involved in a tradeoff relation. Therefore, when 
rate is small, we can both decrease energy and increase rate at the same time. 

• Route with more relays tend to consume less transmission energy. 

• A fundamental limit of the transmission rate exists when relays' number k > 4, due 
to the symmetry of the transmission scheme and the interference from other nodes. 



6.2 Guidelines on the design of routing strategy 

A joint consideration of the transmission rate, path loss exponent, the length of a route 
and the number of relays is necessary, in order to select a path with the best performance in 
terms of EE, BE and latency. The transmission rate and path loss exponent are utilized 
to determine whether TWRC performs better than direct transmission, and whether 
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relay in the middle gives the best performance. If the transmission rate is low, then the 
processing energy is dominant, a route with a small number of relays tends to offer the 
best performance. As rate increases, transmission energy becomes dominant, then a route 
with shorter length and more relays tends to consumes less energy at the cost of higher 
latency. 

7 Conclusion 

In this paper we presented an information theoretical framework to study routing path 
selection for amplify-and-forward two-way relay networks. We formulated bandwidth effi- 
ciency, energy efficiency and latency for routes with different number of relays, assuming 
a simple Hop-by-Hop scheduling scheme. We also determine the optimal power alloca- 
tion scheme that allows a route to consume the minimal energy while still achieve the 
same end-to-end transmission rate. Our theoretical formulations and simulation results 
help provide guidelines towards routing protocol design. Future work of interest includes 
considering routing path selection under the assumptions of multiple frequencies and im- 
perfect scheduling. 
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